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Abstract. A perturbation theory for the Nonlinear Schrodinger Equation 
(NLSE) in ID on a lattice was developed. The small parameter is the strength 
of the nonlinearity. For this purpose secular terms were removed and a prob- 
abilistic bound on small denominators was developed. It was shown that the 
number of terms grows exponentially with the order. The results of the per- 
turbation theory are compared with numerical calculations. An estimate on 
the remainder is obtained and it is demonstrated that the series is asymptotic. 



1. Introduction 

We consider the problem of dynamical localization of waves in a Nonlinear 
Schrodinger Equation (NLSE) [1] with a random potential term on a lattice: 

(1.1) id t ^ = —J [ip(x + l)+tp(x — 1)] + e£V + M 2 1> 

where i\) = \j} (x, t) , x € Z; and {£"} w6 q is a collection of i.i.d. random variables 
chosen from the set fi, with the probability measure /i (e x ) . It will be assumed that 
H{e{x)) is continuous, bounded and of finite support and, additionally, symmetric, 
= M( — £ x)- We assume that exponential localization is known to take place 
for all the energies of the linear problem (when (3 — 0). The decay rate 7 and the 
localization length, £ = I/7, are given for the linear part of the model (1.1) by the 
Thouless formula [ ]. In particular, if fj,(e x ) is a uniform distribution than as was 
found numerically (see Appendix), the function ~f(E) is unimodal. 

The NLSE was derived for a variety of physical systems under some approxima- 
tions. It was derived in classical optics where rp is the electric field by expanding 
the index of refraction in powers of the electric field keeping only the leading non- 
linear term [3]. For Bose-Einstein Condensates (BEC), the NLSE is a mean field 
approximation where the density f3\ip\ 2 approximates the interaction between the 
atoms. In this field the NLSE is known as the Gross-Pitaevskii Equation (GPE) 
[4, 5, 6, 7, 8, 9]. Recently, it was rigorously established, for a large variety of in- 
teractions and of physical conditions, that the NLSE (or the GPE) is exact in the 
thermodynamic limit [10, 11]. Generalized mean field theories, in which several 
mean-fields are used, were recently developed [12, 13]. In the absence of random- 
ness (1.1) is completely integrable. For repulsive nonlinearity (/3 > 0) an initially 
localized wavepacket spreads, while for attractive nonlinearity (j3 < 0) solitons are 
found typically [1]. 

It is well known that in ID in the presence of a random potential and in the 
absence of nonlinearity (/3 = 0) with probability one all the states are exponentially 
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localized [14, 15, 16, 17]. Consequently, diffusion is suppressed and in particular 
a wavepacket that is initially localized will not spread to infinity. This is the 
phenomenon of Anderson localization. In 2D it is known heuristically from the 
scaling theory of localization [18, 16] that all the states are localized, while in 
higher dimensions there is a mobility edge that separates localized and extended 
states. This problem is relevant for experiments in nonlinear optics, for example 
disordered photonic lattices [19], where Anderson localization was found in presence 
of nonlinear effects as well as experiments on BECs in disordered optical lattices 
[20, 21, 22, 23, 24, 25, 26, 27, 28]. The interplay between disorder and nonlinear 
effects leads to new interesting physics [26, 27, 29, 30, 31, 32]. In spite of the 
extensive research, many fundamental problems are still open, and, in particular, 
it is not clear whether in one dimension (ID) Anderson localization can survive the 
effects of nonlinearities. This will be studied here. 

A natural question is whether a wave packet that is initially localized in space 
will indefinitely spread for dynamics controlled by (1.1). A simple argument in- 
dicates that spreading will be suppressed by randomness. If unlimited spreading 
takes place the amplitude of the wave function will decay since the L 2 norm is 
conserved. Consequently, the nonlinear term will become negligible and Anderson 
localization will take place as a result of the randomness. Contrary to this intu- 
ition, based on the smallness of the nonlinear term resulting from the spread of 
the wave function, it is claimed that for the kicked-rotor a nonlinear term leads to 
derealization if it is strong enough [33]. It is also argued that the same mecha- 
nism results in derealization for the model (1.1) with sufficiently large (3, while, 
for weak nonlinearity, localization takes place [33, 34]. Therefore, it is predicted 
in that work that there is a critical value of f3 that separates the occurrence of 
localized and extended states. However, if one applies the arguments of [33, 34] 
to a variant of (1.1), results that contradict numerical solutions are found [35, 36]. 
Recently, it was rigorously shown that the initial wavepacket cannot spread so that 
its amplitude vanishes at infinite time, at least for large enough (3 [37]. It does 
not contradict spreading of a fraction of the wavefunction. Indeed, subdiffusion 
was found in numerical experiments [33, 37, 38]. In different works [38, 39, 40] 
sub-diffusion was reported for all values of (3, but with a different power of the 
time dependence (compared with Ref. [33]). It was also argued that nonlinearity 
may enhance discrete breathers [31, 32]. In conclusion, it is not clear what is the 
long time behavior of a wave packet that is initially localized, if both nonlinearity 
and disorder are present. This is the main motivation for the present work. Since 
heuristic arguments and numerical simulations produce conflicting results, rigorous 
statements are required for further progress. 

More precisely, the question of dynamical localization can be rigorously formu- 
lated as follows: assume the initial state is, ip (x,0) = uo (x) , where uq (x) is an 
eigenstate of the linear part of (1.1) which is localized near x — 0. Then for any 
< 6 < 1, one has to prove that with probability 1 — S (on the space of the 
potentials) 



(1.2) sup e v ^4>(x,t) 

X,t 



< M s < oo 



for some v > 0. 

Rigorous results on dynamical localization for the linear case are well known 
[41, 42, 43]. However, the nonlinear problem turns out to be very difficult to 



PERTURBATION THEORY FOR THE NONLINEAR SCHRODINGER EQUATION. 



3 



handle, even numerically. Consider the case of small (3. There are two possible 
mechanisms for destruction of the localization due to nonlinearity. 

One way of spreading is to spread into many random places with increasing 
number of them. In that case, due to conservation of the normalization of the 
solution, the solution becomes small. But then, the nonlinear term becomes less and 
less important and we expect the linear theory to take over and lead to localization. 
While this argument sounds plausible there is no proof along this lines. 

The second way of spreading is in a few fixed number of spikes that hop randomly 
to infinity. In this case, the nonlinear term is always relevant. It is this (possible) 
process that makes the proof of localization in the nonlinear case so elusive. It also 
precludes a quick numerical analysis of the problem: it may take exponentially long 
time to see the hoping. 

Rigorous results in this direction are of preliminary nature: In [44] it was shown 
that dynamical localization holds for the linear problem perturbed by a periodic 
in time and exponentially localized in space small linear perturbation. In [45] the 
above result was extended to a quasiperiodic in time perturbation. Such perturba- 
tions mimic the nonlinear term: 

2 



(1.3) |V| 2 



Esi 



.1 



where Uj are the eigenfunctions of the linear problem with energies Ej . However in 
other situations time dependent terms may result in derealization [46, 47]. Using 
normal form transformations Wang and Zhang [48] studied the limit of strong dis- 
order and weak nonlinearity, namely, e = J + /3, small. For initial wavefunctions 
with tails of weight 6 starting from point jo, they have proved, that the wave- 
function spreads as following. There exist C — C (A) > and e (A) > and 
K = K (A) > A 2 such that for all t < (8/C) e~ A the weight of the tails of the 
spreaded wavefunction starting from j = jo + K is less that 25. On the basis of 
this result they have conjectured that the spread of the wave function is at most 
logarithmic in t. 

Furthermore, it can be shown that NLSE has stationary solutions 
(1.4) E^ E = {-d xx + e%)^ E + P\^ E \ 2 ^E 

which are exponentially localized for almost all E with a localization length that is 
identical to the one of the linear problem [49, 50, 51, 52, 53]. 

In our previous work [54] we have developed a perturbation theory in (3. By 
considering the first order expansion we have proved that for times of order O (/3~ 2 ) 
the solution of (1.1) remains exponentially localized. A result of similar nature for 
a nonlinear equation of a different structure was obtained in [55]. In the current 
work we consider an expansion of any order, N, in (3. This expansion enables in 
principle the calculation of the solution to any order in N. A bound on the error 
can be computed using only propreties of the linear problem ((3 — 0). Therefore 
this work has the potential to develop into a method for solution of some type of 
nonlinear differential equations. In Section 2 we construct the solution as a series 
in the eigenfunctions of the linear problem. Standard perturbation theory for the 
coefficients does not apply: we encounter small divisor problems and secular terms 
(formally infinite). Removing the secular terms requires the "renormalization" of 
the original linear Hamiltonian by shifting the energies (Section 3). The estimates 
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of the small divisor terms are performed in the spirit of the work of Aizenman- 
Molchanov (A-M) [56]. In Section 4 the entropy problem is resolved by bounding 
an appropriate recursive relation. A general probabilistic bound on the terms of 
the perturbation theory is derived in Section 5 and the quality of the perturbation 
theory is tested in Section 6. In Section 7 the remainder terms are controlled by 
a bootstrap argument. The results are summarized in Section 8 and the open 
problems are listed there. 

In summary, in this work a perturbation theory for (1.1) in powers of (3 was 
developed and bounds on the various terms were obtained. The work is only partly 
rigorous. In some parts it relies on Conjectures that we test numerically. 



2. Organization of the perturbation theory 
Our goal is to analyze the nonlinear Schrodinger equation 

(2.1) id t1 p = H iP + Pm 2 iP 
where Hq is the Anderson Hamiltonian, 

(2.2) H iJj (x) = -J [V (x + 1) + ^ (x - 1)] + e x i) (x) . 

We assume throughout the paper that H satisfies the conditions for localization, 
namely, for almost all the realizations, uj, of the disordered potential, all the eigen- 
states of Ho, u m , are exponentially localized and have an envelope of the form 
of 

(2.3) \u m (x)\<D^ E e^e-^- x ™\, 

where e > 0, x m is the localization center which will be defined at the next sub- 
section, 7 is the inverse of the localization length, £ = and D UiB is a constant 
dependent on e and the realization of the disordered potential [57, 58] (better es- 
timates were proven recently in [59, 60]). It is of importance that -D u , e does not 
depend on the energy of the state. In the present work only realizations ui, where 

(2.4) \D^ e \ < D s>e < 00 

are considered. This is satisfied for a set of a measure 1 — 6, since (2.3) is false 
only for a measure zero of potentials. 

2.1. Assignment of eigenfunctions to sites. It is tempting to assign eigenfunc- 
tions to sites by their maxima, namely, is the eigenfunction with energy Ei and 
a maximum at site i. This assignment is very unstable with respect to the change 
of realizations. This is due to the fact that the point where the maximum is found, 
which is sometimes called the localization center, may change as a result of a very 
small change in the on-site energies {e x }. To avoid this, the assignment is defined 
as the center of mass [61], 

(2.5) x E = y~] x \u E , (x)\ 2 

Definition 1. The state u Ei is assigned to site i if i — [xe] ■ If several states are 
assigned to the same site we order them by energy. 
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2.2. The perturbation expansion. The wavefunction can be expanded using the 
eigenstates of H as 



(2.6) 



i/j (x, t) = ^2 c m (t) e lEmt u m {x). 



For the nonlinear equation the dependence of the expansion coefficients, c m (t) , is 
found by inserting this expansion into (2.1), resulting in 



(2.7) id t c m e- iEmt u m (as) = H ^ c„ 



e- iE ^u m (x) 



■0 



E 



1 ' nit- 



-iE m t 



^ ] c m 3 e lEm3t Umz ( x ) ■ 



Multiplying by u n (x) and integrating gives 
(2.8) id t c n = P J2 V n Hl 



m 2 m 3 * i(E mi +E n ~E m2 -E m )t 



mi ,ni2 ,m3 

where F™ 1 " 12 " 1 ' is an overlap sum 

ymim 2 m 3 



(2.9) 



= ^ u n (%) U mi (x) U m2 (x) U m3 (x) . 



By definition V^ lim2m3 [ s symmetric with respect to an interchange of any two in- 
dices. Additionally, since the u n (x) are exponentially localized around x n , V™ 1 ™ 2 " 13 
is not negligible only when the interval, 

(2.10) 5m = max [x n , x mi ] - min [x n , x mi ] , 

is of the order of the localization length, around x n , 

(2.11) 

|^mim 2 m3| < £)j e e(\x„\+\x mi | + |s m2 1 + 1^3 |) _ e —y(\x-x n \ + \x-x mi \ + \x-x m2 \ + \x-x m3 

x 

< Dj £ e £ (l a: "l + l a: '"i \ + \ x ™2 \ + \ x m 3 |) e - <T 3 - 1 (\x n -x mi \ + \x n -x m2 | + |x„-a m3 |) x 
x , e -e'(\x-xn\+\x-x mi \+\x-x m2 \+\x-x m3 \) 

X 

< Vg' e ' e E ^ x '^ + \ Xmi \ + \ x ™2\ + \ x ™ 3 \) e -H~'- e ')(\ x <^ x ™i \+\ x "- Xm 2\+\ Xn - Xm a\) 
Here we have used the triangle inequality 

(2.12) (\x - x n \ + \x- x mi \) + (\x - x n \ + \x- x m2 \) + 

~t~ y\x x n | -(- \x Xjyi^ |) ^ \x n x mi | -I- \x n x rri2 1 ~t- \x n x m3 1 

to obtain the second line. Our objective is to develop a perturbation expansion of 
the c m (t) in powers of (3 and to calculate them order by order in (3. The required 
expansion is 

(2.13) c n (t) = C W + p c W + +... + p^d^ + P N Qn, 

where the expansion is till order (N — 1) and Q n is the remainder term. We will 
assume the initial condition 



(2.14) 



c„ (t = 0) = S n0 . 
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The equations for the two leading orders are presented in what follows. The leading 
order is 

(2.15) 4°) = 5 n0 . 

The equation for the first order is 
(2.16) 

.■/J, r (l) - \" V"»im 2 m 3 *(0) (0) (0) i(E n +E mi -E m2 -E m3 )t _ t^OOO i{E n -E a )t 
mi,m2,m3 

and its solution is 



(2.17) 



(1) = yOOO 



1 _ e i(E n -E )t 



E n — E 

The resulting equation for the second order is 

1 m 2 m 3 

v n 

mi,m2,m3 



(2.18) id t c^ = J2 C im2m3 cj 1 1) ^^e i(£ " +£ '" 1 ^ m2 ^ m3)i + 

4-2 V y™i™2m 3 „*(0) (1) (0) i(E n +E mi -E m2 -E m3 )t 
' " / j v n "uii u m 2 "13 



mi ,m 2 ,m3 

Substitution of the lower orders yields 

, 1 _ p -i(E m -E a )t 



(2.i9) zd t c! 2 > = ]Tv;r 00 v; 

m 

' 1 — e i ( Em ~ E ° >t 



+2 



E 



-Em — E 
ymOOT/000 



-Em — E a 
,i(E n -E m )t 



e i(E n +E m -2E )t + 



,n E m-E 

+2 (^KE n -E m )t - e »(^n-s )t^ 

ymOOyOOO 



^ e i(E n +E m -2E a )t _ e i{E n -E )t^ _|_ 



E 



're ' m 
E m — E 



e i(E n + E m -2E )t _ £ e i(E n -E )t _|_ 2e *(-E„-B m )t 



We notice that divergence of this expansion for any value of /3 may result from three 
major problems: the secular terms problem, the entropy problem (i.e., factorial 
proliferation of terms), and the small denominators problem. 

3. Elimination of secular terms 

We first show how to derive the equations for c„ (t) where the secular terms are 
eliminated. 

Proposition 2. To each order in (3, ip (x, t) can be expanded as 

(3.1) i>{x,t) = J2 c n{t)e- iE '» t u n {x) 

n 

with 

(3.2) E' n ^E^ + ^+^E^ + --- 

and En^ 1 are the eigenvalues of H , in such a way that there are no secular terms 
to any given order. The E' n are called the renormalized energies. 
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Here we first develop the general scheme for the elimination of the secular terms 
and then demonstrate the construction of E' n when the c„ (t) are calculated to the 
second order in (3 (see 3.19,3.18). 

Inserting the expansion into (2.1) yields 

(3.3) 

*E idtcm - iE' mCm ] e- lE '^u m (x) = ^E^e"^'^ (x) + 

m 77i 

+/3 E C*n 1 c m 2 C m3 e( Em i~ Em 2~ Em ^ t U„ ll (x) U m2 (x) U m3 (x) . 

771177127713 

Multiplication by u n (x) and integration gives 
(3.4) 

where the ^"^^^ are given by (2.9). Following (2.13) we expand c„ in orders of 
(3, namely, 

(3.5) c n = ci°) + /?c« + (3 2 c^ + •••,. 

Inserting this expansion into (3.4) and comparing the powers of (3 without expand- 
ing the exponent, produces the following equation for the k — th order 

(3.6) 



E v ™ 17 



k—1 k — 1— r k—l—r—s 

E E E 

|_r=0 s=0 ;=o 



Jr)*Js) (I) 

7711 7712 7713 



Note that the exponent is of order O (1) in /3, and therefore we may choose not to 
expand it in powers of (3. However, it generates an expansion where both Em and 
ci^ depend on (3. For the expansion (2.13) to be valid, both E m and c„ should be 
O (1) in (3, this is satisfied since the RHS of (3.6) contains only cL^ such that r < k. 
Namely, this equation gives each order in terms of the lower ones, with the initial 
condition of cl°' (t) = S n Q . Solution of k equations (3.6) gives the solution of the 
differential equation (3.4) to order k. Since, the exponent in (3.6) is of order O (1) 
in (3 we can select its argument to be of any order in (3. However, for the removal 
of the secular terms, as will be explained bellow, it is instructive to set the order 
of the argument to be k — 1, as the higher orders were not calculated at this stage. 
Secular terms are created when there are time independent terms in the RHS of 
the equation above. We eliminate those terms by using the first two terms in the 
first summation on the RHS. We make use of the fact that Cn^ — <5 n o an d Cn can 
be easily determined (see (3.9,3.12)), and used to calculate E^2o an d that 

(k) 

eliminate the secular terms in the equation for c„ , that is 

T/000 

(3.7) E^c^+E^c^^E^Sno + E^a-Sno) " 



K-E' ' 
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where only the time-independent part of c„ ' was used. In other words, we choose 
E^ and £fc so that the time-independent terms on the RHS of (3.6) are elim- 
inated. Eq^ will eliminate all secular terms with n = 0, and E^ ^ will eliminate 
all secular terms with n ^ 0. In the following, we will demonstrate this procedure 
for the first two orders, and calculate Cn \ and obtain an equation for d 2 '. 
In the first order of the expansion in j3 we obtain 

(3.8) ift<£> = + C im2m3 ^)^cJS) e '(«r^-^)' 

= -EW5 n0 + V° 00 e i «- E °>. 
For n = the equation produces a secular term 

(3.9) tdtc^ = -E™ + V Q 000 

Setting 

(3.10) = V 00Q 
will eliminate this secular term and gives 

(3.11) c« = 

For n^O there are no secular terms in this order, therefore finally 



(3-12) c« = (1 - 6 n0 ) V° 



K - K 



where to this order E' n = E n and E' = Eq- 

In the second order of the expansion in (3 we have 

(3.13) 

mim2m3 

= -E^6 n0 ^ + E C 1100 (c^V^^" 2 ^ + 2c«e«-< 1 > 

mi 

For n = it takes the form 

id tC M = -4 2) + E K. m0 ° (c^V^"^)* + 2cWe'( £ i-^) f ) . 
m 

Substitution of (3.9) and (3.12) yields 

T/mOOT/000 
„(2) _ „(2) , 10 V _r 



(3.14) id t c^ = + V K ° "\ [f 1 - e-'C^-^)*) e <^-^)*+ 

+2 (l - e't^™-^)*) e< E °- E '™y 



T/mOO-t/0 00 / 

£ (2) + V- Vv_Kn_ Li{E' m -E' Q )t + 2e i{E',-E' m )t _ g \ ^ 

, - E' E n V / 
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and the secular term could be removed by setting 

(3.15) ^ = -*E\^r- 

For n/Owe have 



(3.i6) id tC ™ = -e^v; 







m 



E&V 



i ii id 



ymOOyOOO 



We notice that the second term in the sum produces secular terms for m = n. Those 
terms could be removed by setting 

pMl/OOO 9T/n00T/000 

(3.17) - En Vn , + 2Vn \ =0 n + 

E 'n ~ E E n~ E 

E{ 1) = 2 yn00 n _^ 

To conclude, up to the second order in f3 , the perturbed energies, which are required 
to remove the secular terms, are given by 



(V°°°\ 

(3.18) E' n = + PV^ (2 - 6 n0 ) - 3[3 2 5 n0 ^ ™ 



E ™ ~ E 



1 > 
o 



and the corresponding correction to is 
(3.19) 



n = 
n ^ 



Note that in the calculation of c„ to higher orders in 0, a secular term of the 
order /3 2 will be generated for n ^ . Secular terms with increasing complexity 
are generated in the cancellation of higher orders, however, as demonstrated by 
(3.7), secular terms are removed with the same <% and Cn which are presented in 
(2.15,3.12). 

In the next section, the entropy problem will be studied. It will be shown that 
the proliferation of terms in the expansion is at most exponential. 
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4. The entropy problem 



Since the time dependence of all orders is bounded (excluding the secular terms) 
we can bound each order of the expansion by 



(4.1) 



-,(0) 



-.(2) 



— $n0 



(li)U 



V, 



E 



1-e 



< 2 



V 



000 



dt 



/T/m00T/000 \ . . 
, / v n V_m \ f e i{E' n +E' m -2E' )t' _ ^{E^-E'^t' + 2e K E n- E m ) t ) 

V E 'm ~ K J ^ ' 



It/toOOI |t/000 I 
o \ " \_n II ™ 
- ' | R/ _ R'| 



et cetera. However, for convergence for a finite but possibly small /?, it is essential 
that the number of terms on the RHS of (4.1) will not increase faster than expo- 
nentially in k, e.g. not as kl, where k is the expansion order. Next we will show 
that the number of terms indeed increases at most exponentially in k. 

We will designate the number of different products of order k oi V's by R\. (on 
top of it there is still a number of non vanishing terms in the sums over m, that 
will be estimated in the next section). By replacing each c^} in (3.6) by Ri (the 
integration with respect to time multiples the number of terms by a factor of 2, cf. 
(4.1)) we deduce a recursive expression for Rk 



(4.2) 



Rk — 



k-1 k-1-1 

E E 

r=0 1=0 



R r RlRk-l-r- 



R\ = 1 Rq — 1. 



In order to find an upper bound on Rk we examine the structure of the products 
of V's we notice that each product could be uniquely labeled by a vector of zeros 
and m^s 
(4.3) 



V" 



m, 1 00 ^/ r ' m 2 m 3 



■ ■v u 



iOt/000 



{mi,0,0,m 2 ,m3,0,- • • , 0, m fe _i, 0, 0, 0, 0} , 



where the number of different summation indices rrii is (fc — 1) and the length of the 
labeling vector is 3fc. Since in each vector the last three elements should always be 
zeros the number of different configurations of this product is the number of ways 
to distribute (k — 1) m's in 3 (k — 1) cells (superscripts), namely, R^ < ( 3 ^Zp)- 
This is only an upper bound, since there may be some additional constraints, for 
example, the first three elements in the vector should never be all zeros. Subtracting 
the cases when all three first elements are zero ( 3 [ k ~P) we obtain the bound 



(4.4) 



Rk < 



3 (k - 1) 
k-1 



(TP) 

3(fc-2) 
k-1 



This bound has the following asymptotics 



(4.5) 



lim — In 

fc-+oo k 



3fc 

k - 



3fc 
k- 



i 27 

ln T' 
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namely, 

(4.6) R k < e fcln ¥ < e 2k . 

From Fig. 4.1 we conclude that this bound is a tight bound of the exact solution 
of the recurrence relation. This bound shows that the number of terms in the 




Figure 4.1. The solid line denotes the exact numerical solution 
of the recurrence relation for Rk and dashed line is the asymptotic 
upper bound on this solution (e 2k ) . 



expansion increases at most exponentially in k and therefore there is no entropy 
problem. 



5. Bounding the general term 



As clear from (3.6) after the subtraction of all the secular terms in the preceding 
orders the differential equation for the fc-th order term is 

(5.1) 



fe-i 



«). 



s=0 



E k 



mim2m.3 



k— 1 k—l—r k—l — r—s 

E E E 



Jr)*Js) (l) 

Kj jTl± ^"7712 TT13 



r=0 s=0 



1=0 



KK+E' mi -E' m2 -E' m3 )t 



with the initial condition of c„ = <5„o and the first term on the RHS is designed (see 
Section 3) to eliminate all the time-independent part of of RHS of (5.1). Following 
the construction of the lower order terms in the preceding section by a repeated 
application of (5.1) the structure of the general term in the expression for a given 
order k can be obtained. Note that the structure is similar to the structure of 



Cn 



(I) 



The main blocks of the structure take the form 



(5.2) 



C 



mi m2fTi3 



v. 



m 117121713 



where {E'} m denotes a sum of eigenenergies (shifted so that the secular terms are 
removed, see (3.2)) that may depend on the summation indices mj. Then any term 
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of order A: is a product of k factors of the form (5.2) and k — 1 summations over 
the indices m, 



k terms 



/r o\ >m 1 m2m 3 >m 4 m5m e _ _ _ >000 ^000 

and following the last section there is an exponentially increasing (in k) number of 
such terms. In order to bound the general term of order k we will first bound one 
typical block, namely, 

ym 1 m 2 m 3 

( ' ^ ~ E' n - {E'} mi 

where {E'} m is some sum of Ej. To bound (5.4) we will bound separately the 
denominator and the numerator. Using Cauchy-Schwarz inequality, 

/ \ 1 I 2 

(5.5) <IC im2m T> < ( -, * - , 2a > (|yr TO2ro3 l 2s ' 

\\E' n -{E'} m ,\ I X 
where < s < |. 

Conjecture 3. For the Anderson model, which is given by the linear part of (1.1), 
the joint distribution of R eigenenergies is bounded, 

p(E 1 ,E 2 ,...,E R )<D R 

where Dr oc R\ < oo. 

The conjecture is inspired by Theorem (3.1) of the recent paper by Aizenman and 
Warzel [62], If one assumes that with probability one the profiles of the eigenfunc- 
tions, namely, the squares of the eigenfunctions, which correspond to the eigenener- 
gies {Ei}f =1 are substantially different such that a (as defined in Theorem (3.1) of 
[62]) is bounded away from zero, than taking the intervals Ij = dEj one finds that 
the joint probability density can be bounded by Dr oc It is not known how to 
prove that for the Anderson model the profiles of the wave functions are distnict 
and how to quantify this. However, it is reasonable to assume distinctness since 
different eigenfunctions are localized in different regions and therefore are affected 
by different potentials. There are double humped states (consisting of nearly sym- 
metric and antisymmetric combinations of two humps), which have approximately 
the same squares, and therefore are natural candidates for states that may result in 
violation of the conjecture. Nevertheless, those states are very rare and the differ- 
ence between their squares is exponential in the distance between the humps. For 
this it is crucial that many sites are invloved (therefore the counter example (2.1) 
of [62] is not generic). If the energies are assigned to specific locations than the 
factorial term could be dropped, namely, Dr oc a~ R . This is due to the fact that 
specific assignment of energies chooses one of the Rl permutations, mentioned in 
[62]. 

R 

Corollary 4. Given < s < 1, for f = ^^c^E ik , where Ck are integers (and the 

k=l 

assignment of eigenfunctions to sites is given by Definition 1) the following mean 
is bounded from above 

(5-6) (~}<D R <oc. 
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where D R oc Dr. 
Proof. By Conjecture 3, 



(5.7) 



(i/r 



P (E 1 ,E2,...,E R )dE 1 dE 2 ---dE R 



R 



k=l 



< Di 



dE 1 dE 2 ---dE R 



R 

E 

k=l 



CkEi, 



changing the variables to {/, E 2 , E3, . . . , Er} gives 



(5. 



(i/r 



< D, 



dE 9 



■dEj. 



— A 



(E) l/l 



where |ci| is the Jacobian and 2A is the support of the energies. Due to the fact 
that f (Ei) is linear the multiplicity is one. Since the integrand is positive we can 
only increase the integral by increasing the domain of integration of /. Designating 
by /oo the maximal value of /, 



(5.9) 



1 

17P 



< 2D R A 



-1 fa 



Dj 



□ 



Conjecture 5. In the limit of R — > 00, for < s < 1 and for f 



R 



Cfe are integers ( and the assignment of eigenfunctions to sites is 
1) 

.i/r/ r s/2 



/c=i 



c k E ik , where 
by Definition 



(5.10) 



For large R the sum, / 



R 

E 

k=l 



CkE ik , can be effectively separated into groups of 



terms that depend on different diagonal energies, Ej. Therefore by the central limit 
theorem, / is effectively a Gaussian variable with (/) = and ^/ 2 ) = a 2 R, where 
a 1 is some constant. Therefore, 
(5.11) 



1 



2 

V2na 2 R 



df _ 
f s 



f 2 /(2<r 2 R) _ 



df _ 
f s 



/ 2 /2 



Conjecture 3 and Corollary 4 were tested numerically for lattice size 128, s 
and the uniform distribution 



(5.12) 



M Q 2A 



\s(x)\<A 
\e(x)\>A' 



with A = 1. The results are presented in Fig. 5.1 for R < 10. For R < 3 all 
the combinations of the energies were used and the result is an average over all 
the combinations. For R > 4 only a partial set of combinations of cardinality 10 4 , 
chosen at random was used. For large R the decay is as R~ s / 2 in agreement with 
Conjecture 5. The above calculation was repeated for the case where the Ei are 
replaced by the renormalized energies E[. The calculation can be performed only 
to the order /3 2 with the help of (3.18). The results are also presented in Fig. 5.1 
for (3 = 0.1 and (3 = 1, 
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0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2. 

In(R) 



Figure 5.1. The logarithm of '' 2 j as a function of the 

logarithm of R. The lines designate denominators with (3 = 0, 
with the solid line (blue) is for J = 0.25 the dashed line (green) is 
for J = 0.5 and the dot-dashed line (red) is for J — 1. The solid 
circles and the squares are data with (3 = 1, and E' n calculated up 
to the second in /3, such that different colors represent different J, 
in the similar manner as for the lines. The solid squares are for 
parameters similar to the ones with the solid circles, but with the 
restriction that at least one of the states that corresponds to E' n 
which is localized near the origin. 



Conjecture 6. Corollary 4 and Conjecture 5 hold also if the Ei are replaced by 
the renormalized energies E[. 

The reason is that the various renormalized energies are dominated by different 
independent random variables The numerical calculations support this point 
of view. In what follows Corollary 4 and Conjecture 6 (and not Conjecture 3) are 
used, and these were tested numerically (Fig. 5.1). 

Using the bound on the overlap sum (2.11) and Corollary 4, 
(5.13) 



ye,s 



s(|a„|+|a: mi \ + \x m2 | + |z m3 |) - i [ 1 -e')s{\x n -x mi \+\x n -x m2 \+\x n - 



(\Cn im2m3 \ S )s< D S 

this proves the proposition: 

Proposition 7. For some S,e,e' > and < s < |, 

(5.i4) (ic nim2m3 r>5 < ^v^'^'+^i^i^M'^o^.k- 3 



where F' = Dx 
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Using the Chebyshev inequality, 
(5.15) Pr(|x| > a) < (\x\) /a, 

where x is a random variable and a is a constant, one finds 

Corollary 8. 

(5.16) 

p r ^mim 2 m 3 | ■> pll/s^ \x n \ x mi \ ) £ - § (is' -77) Y.i\xn-x m ^ < 



where F' = Dx 



A general term in the expression for different orders of d^' is given by the form of 
|^ro.im 2 m 3 | j^msme | . . . £000 ^ ( ie>; it conta ins (k - 1) summations over indices 

which run over all the lattice. First we construct a general procedure to bound a 
product of k, ( 's. A product of two £ 's is bounded by 



) <(£ic 

5 \ mi 



where (-} s denotes an average over realizations where (2.4) is satisfied. 
Using the Cauchy-Shwarz inequality 

1 /2 



1/2 
S 



setting < s < 4 (notice, that s < | and not s < |, due to (5.5)) and inserting 
the bound on the average ^|£™i m 2 m 3 | 2s ^} from Proposition 7 gives 



(5.17) 



E„m 3 /-Tl4? 



< F' exp 



Ei 



e 3 



"(I 



> e 2es l^ilexp-s^^ 



where we have used the inequality 



En) <£n 

V i / i 



< s < 1. 



Using the triangle inequality in the same manner as in (2.12) 

6 6 6 

(o.l8) |x n X-mi | ~t~ ^ ^ |^mi ^m; | 77 ^ ] | 77 ^ ] I^JTii 



i=4 



i=4 



i=4 
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we get 



(5.19) 

/ V ^rtnimjffls /-mil 



() 



<F'exp £s [\x n \ + J2\x mt \j e-^'CK-^l+k"--!): 
7-e' 



x exp 



= F"exp 



9 



5^ |X„ 



2 = 4 



E 



c 2BS F"".,e 



S£ | |x n | H - ^ ^ I ^T^Xi I J q ^ ^mj 



where F" {^,e f ,s,e) = F 



mi 



-2^s £ |x mi -x m< | 



< oo, in the following 



also other convergent sums of this type will be denoted by F". 

If the term we consider is a term in the perturbation expansion it should include 
some factors £™i m 2«i3 with some m, = 0. A simple example is where x mi = x m2 = x m 



(5.20) 



(|CT) <F'e-<i- e '- £ >^, 



where the bound was calculated using Proposition 7. The product should terminate 
with a term of the form therefore a term like (5.17) is a part of a product of 
the form, 



(5-21) Yl 



\^m 1 m 2 m^\ \/-m 4 m, 5 m G \ I ^000 I I aOOO I I /-000 I I a000 I I ^000 I 

I | ami I I ^7712 | | =>"T<3 I I =>7Tl4 | | ^5 | 1^7716 | 



To bound it we use the generalized Holder inequality, 



l/k 



Applying it yields, 



tr nn\ \ " I l/-m 1 m 2 m 3 \s \^m4,m s m 6 \ s I a000 I s I a000 I s I a000 I s I ^000 I s I a000 I s v 

i y / ^ ^l^n I l^mi I |^m2 I I ^7713 | |^m4 | | | l^mg I j 

/ 6 \ V7 



< E (ic ir 



1 7 s 



n ic 



i=2 



^000 



7 s 
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aOOO 




Figure 5.2. An example of a graph that is used to construct 
the general term. The graph describes an 8-th order term, 

Smi ^7773 Sm2 ^7774 >mg We W7 • 



setting < s < jj and inserting the bounds on the averages (5.17) and (5.20) gives 
(5.23) 



^000 

■57716 



\ Ami m 2 7713 Am 4 m5m6 aOOO aOOO aOOO aOOO 

/ ^ Ss77 ^7772 ^ 7773 ^7774 ^7775 ^77 

<F'J2 expes ( \x n \ + £ \x mi \ J e - 3t ¥ i -(l»»-*» a |+K-*»8|)x 

{777,} \ j=2 / 

, 6 6 

7-e 



x exp 



-^y^ \x n - x mi \ - (7 - e' - e) |a 

7=4 7 = 2 

, 6 6 

1 q i x " ~ ^1 - (7 - £ ' - 2e ) 1 



x exp 



7=4 7=2 
where {mj stands for a sum over all the m,. Using the inequality, 

6 

(5-24) £(| l)>3| ^77 1 1 



7=4 



we get 
(5.25) 



V" ^777 1 777 2 7773 ^777477757776^000^000^000^000 /■ 
/ j ^77 ^ 7771 ^7772 ^ 7773 ^7774 ^7775 V, 

A m i} 



000 

7776 



< F > e -b-e'-e)s\ Xn \ J- eX p-2 S l 1 -^ - s)j2 1^1 

{777,} ^ ' «=2 



^' e -(7-e'-e)sk»|_ 



Let us study the form of a general graph (c.f. Fig. 5.2). 
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It can be described as a tree starting from n, the "root" and four types of branch- 
ing points, where q (0 < q < 3) branches continue while 3 — q branches terminate. A 
branch which continues is associated with a value m, ^ 0, while a branch terminates 
with mi — 0. In the above, bounds on branches with q = and 3 are calculated, 
and are given by (5.14) and (5.20), respectively. The bounds for q = 1,2 follow 
similarly from Proposition 7. Along each "bond" from the "root" to the "leaves" a 
term (™ im!m3 is multiplied. At a point from where q branches continue the ex- 
ponent is reduced by a factor of q. In other words, with x n and x m are connected 
by a path that is crossing I branching points with ratios qi , . . . , qi the bound on 



the product of the zetas contains a factor exp - 



i K. 



where q = The 



total number of branch ends is also q. To terminate all the branches q factors of 
the form (bounded by (5.20)) should be multiplied resulting in the a term that 
multiplies the bound on a sum of the form 



(5.26) 



E exp 

{m,} 



< 



exp - 



(7-e') _ 

q 

{l-e 1 ) 

q 

-sq\x n \ 

{mi} 



Ei 



s E \ x * 



exp- 1 - - (-f-e')s 



exp - 



( 7 -e')sJ2\xm i \ = S s 



-(-y-e')s\x n \ 



restoring the original convergence rate. As the product consists of k terms the 



evaluation of 



sk- 



is required for the use of the Holder inequality. Therefore 



it is required that < s < l/2k. 

Lemma 9. For a given k (the number of <^'s), S,e,e',f]' > and < s < A- 



(5.27) 



Ec 1 



•C 



000 

mjv-i 



or using the Chebyshev inequality (5.15) 
(5.28) 



Pr 



morris /-m^m^iriQ 



•C 



000 
mjv-i 



> 



1/., 



_^ 7 _ e _ e '_^Q| Xn | 



< e 



-77 s\x n I 



where F$ is a constant which is built iteratively by the construction demonstrated 
in (5.23) and (5.25) and is proportional to Dg. 

It is of importance, that any product with the same number of zetas has the 
same bound with the same probability. This allows us to bound Cn by counting 
the number of different configurations, Rk, of the product for a given k and then 
multiplying it by the bound of each product. This proves the theorem: 



Theorem 10. For a 



k and 5, e, e', 77' > 



(5.29) 



Pr 



> 



ck 2 +c'k -(7— e-e'-)7')|x„| 



< e" 
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where Fi k ^ which is proportional to Dg and c and c' are constants. 



Proof. Using Lemma 9 and summing over configurations denoted by i c one obtains 

the bound 

(5.30) 

«* r 

(fe) \ _ / \ ) \ /■m 1 m2m 3 /■m i m 5 m & _ _ _ ^000 

/ v ) / y Wi 
ic=l I W 




i c =l 



E/ fc mim2m3/ fc m4m5m6 > > _ ^000 




< R k F^ k) e-^- e - e >^. 



Using (4.6) 

(5.31) 
or 

(5.32) Pr 



,(fe) 



< e 2fc j pC fc ) e -(T-£-e') s l ;c "l 



> e C ' /5 ' ( F (k) ) 1/S e -h-c-t'-v')\x n \ e 2k/s\ < e -c' e - V 's\x n \ 



Choosing the largest s produces a bound where c, d are constants and is a 
constant which is built iteratively by the construction demonstrated in (5.23) and 
(5.25) and is proportional to Dg. □ 

Remark 11. From (4.6) one sees that c ±2 2 and later we set c' = cN . 

6. Numerical results 

In this section we will check how well the perturbation series up to the second 
order in j3 approximates the numerical solution of (1.1). For this purpose we use the 



expressions for c4 and Cn' which were obtained in (3.12) and (3.19), respectively, 
and also the expression for the renormalized energies E' n up to second order in (3 
which are given by (3.18). We use the perturbation expansion up to the second 
order in j3, namely 



J 2 ') 



(6.1) 



where we took 



_ JO) 



(2) 



/3=0 ' 



„(2) 



(2) 

in order to keep only the contribution of c„ up to the order 



1 . To compare, we plot the real and the imaginary parts of both the numerical 
solution of (1.1) with the distribution (5.12) and the perturbative approximation 
c„. From figures 6.1 and 6.2 we see that the correspondence between the numerical 
solution of (1.1) and the perturbative approximation is good for times < 50 for 
(3 = 0.1 and times < 200 for j3 = 0.01. Additionally, the correspondence of the 
central site, cq, which is used as the initial condition, c n (t = 0) = S n o, is much better 
than the correspondence of the neighboring sites. A possible explanation for this 
could be that the nonlinear perturbation is more pronounced at the states n with 
c n (t = 0) = 0, this is due to the fact that for (3=0 those sites are unpopulated (zero) 
for all times, resulting in lower signal to noise ratio. To examine the convergence 
in time we will define a time, t*. For times t < t* the relative difference between 
the L 2 norms of the exact and the approximate solutions, d 2 (t), defined bellow, is 
less than 10%. It is instructive to introduce the following definitions. 
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= 0.1 site = 




Figure 6.1. In the first and the second rows of the figure are pre- 
sented the imaginary and real parts, respectively, of the numerical 
solution of (1.1) (blue, solid) and the perturbative approximation 
(red, dashed) as function of time for Co (left column) and ci (right 
column). In the third row the relative difference, di (t) , as de- 
fined in Definition 12 is plotted. The parameters of the plot are, 
J = 0.25, A = 1 and (3 = 0.1 and lattice size of 128. See (1.1) for 
the definition of the constants. 



Definition 12. The relative difference between the exact and the perturbative 
solution at site n is defined as 

2\c n (t)-c n (t)\ 

lU K(t)\+\c n (t)\- 



Definition 13. t* is a time until which the relative difference between the exact 
and the perturbative solution, d,2 (£), is less than 10%. 



(6.2) da (t) = En K W < 01 

E»K (*)l 

In figure 6.3 we see that as (3 becomes smaller the time for which the expansion to 
second order in (3 is close to the exact solution (within 10%), t* , increases. 
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Figure 6.2. Same as Figure 6.1 but for (3 = 0.01. 




Figure 6.3. t* (see text) as a function of the inverse nonlinearity 
strength, /3 _1 . 

7. Bounding the remainder 

In order to bound the solution we have to bound, Q n , the remainder of the 
expansion (2.13) 

(7.1) c n (t) = 4°) + /3c« + ^ + ■■■ + p N -'cl N -V + (3 N Q n . 

This is achieved applying the bootstrap argument to the remainder. Substituting 
the expansion (2.13) into the equation (3.4) and writing a differential equation for 
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the remainder gives 
(7.2) 



N-l 



id t Q n = Yl E /3 fc+ ' +9+1 - N K" im2m3 ^ ) c«cWe<< +£ ™ 1 -^ 

m 1 m 2 m 3 N-l<k+l+q<3(N-l) 
N-l 

+ E E ^r^Me'^-'-^- 1 ^)' 

miin 2 m 3 0<fc+(<2(JV-l) 
N-l 

+ 2 E E /3 Z+fe+ V™ im2m3 m 3cW*cWe< £ " + ^ 1 -^ 2 -<3)* 

roim 2 m 3 0<k+l<2(N-l) 

+ 2 ^ ^ l ^ +fe+ V-^^3Q* ii Q m2C (fe) e <K + < 1 -< 2 -< 3 )* 
TOim 2 m 3 fc=0 

+ E E^ +fc+1 ^ im2ms ^^^^3 eiK+<i " B;,a " B;,8)t 

mim2m3 fe=0 



where the sums over orders are understood as follows. By Sjv-i<fc+z+g<3(N-i) we 
mean ^2kTq=o with the constraint N — I < k + I + q . Integrating and using the 
fact, Q n (i = 0) = 0, 



(7.3) 

\Qn\ 



N-l 

< V V pk+l+q+l-N \yrn im2 m 3 \ gup 

* r-^ o<t'<t 





r *(fe) 




C<'> 




r (9) 








L m 2 




°m 3 



mim 2 m 3 N-l<k+l+q<3(N-l) 
N-l 

+ E E /3 ;+fc+i ic 11 " 

roim 2 m 3 0<fc+i<2(N-l) 
N-l 

_|_ 2 ^ ' ^ ' pl+k+l ^ym 1 m 2 m 3 

m 1 m 2 m 3 0<k+l<2(N- 1) 
N-l 

+ 2 ^ ^iV+fc+l^™^! gup \\Q* \\Q 

0<t'<t L 



I sup 

0<t'<* 



I sup 

0<t'<t 



Qmi I 
|Qm 3 



mim 2 m 3 k=0 
N-l 



+ E V/? Ar+fc+1 |K" im2m3 l sup |Q TO2 ||Q 

0<t'<t 1 

mim 2 m 3 fe=0 
_|_ ^2N+1 y ' |-^mim 2 m 3 



,(fc) 



,(fe)» 



mim2m3 



SUp [|Q^J |Qm 2 l |Q 
0<t'<t 



where i is the upper limit of the integration. 

Since Q n (t) is continuous and Q n (t = 0) = 0, for small i, Q„ (i) <~ t ■ S\. 
Therefore we can always find a sufficiently small r > 0, such that 



(7.4) 



\Q n (t)\<2t-S 1 , 
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where 
(7.5) 



N-l 



m 1 m 2 m 3 N-l<k+l+q<3(N-l) 

Assume, that there is some time, t, for which, 



sup 

o<t'<t 



(7.6) \Q n (t)\ >2t-S u 

than since Q n (t) is continuous and (7.4) holds there is a time, to, where 

(7.7) |Q«(to)| = 2t -Si. 

Inserting this equality into the inequality (7.3), we get an interval < t < to for 
which (7.4) holds (see (7.25)). We proceed by bounding Si and other sums of (7.3). 

Using the Theorem 10 obtained in the end of the Section 5, we can bound, Si. 
The inequality is violated with a probability found from (5.29). We start with the 
bound 



JV-l 



(7. 



si< E E 



k+l+q+l—N Iymim2m3 



{mi} jV-l<fc+Z+g<3(iV-l) 
xe c'(k+l+q)+c(k 2 +l 2 W) e -(-l-e'- V ')(\x mi \ + \x m2 \ + \x m3 \) 

we find from (2.11) 
(7.9) 



^ |y^l™2™3| e -(l-£'-v')(\x mi | + |^m 2 | + |^m 3 | 

mim2m3 



x e 



mim2m3 

-\(l-e')(\x n -x mi \+\x n -x m2 \+\x n -x m3 |) ( 7 - 7) '- e ')(| Xmi | + | Xro2 | + |x m3 |) 



mi 7712^3 

for 2 < (7 — rf) . Therefore for rf sufficiently small, substituting back we get 
(7.10) 

N-l 

Si < C , J' E,e "'' 7 'e~( 7 ~ e_e ')l X ™l pk+l+q+l-N e c(k 2 +l 2 +q 2 ) e c' (k+l+q) 

N-l<k+l+q<3(N-l) 
= C S (N) e 6ciV 2 e -(7-e-s')k„l ; 

where we used d = cN, and Cg (N) = O (iV 3 ). For the bound (7.10) to be violated 
at least one of the Cm' has to satisfy the inequality (5.29) and the probability for 
this is bounded by e~ c e~ n \ x ™\/ k . Therefore the probability that (7.10) will be 
violated is bounded by 
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c (l) 




M 




Cm i 




L m2 




c m3 





where we have expanded the exponent e~ v ^ k using the fact that for large N the 
sum is dominated by terms with k ^> 1 and rj /k <C 1. Setting d = N provides the 
convergence of the probability with the expansion order, i.e., 
(7.12) 

(V 1 sr--N-l ak+l+q+l-N It/m-i m 2 m 3 1 

> C s (N) ,ficN» e -(~t-*-*-J)\*~\ 

Now we turn to find the point to defined in (7.7). To bound other expressions in 
(7.3), we use 

(7.13) |Qm(to)| = 2t ■ Si < M{t a )e- V ^\ 
where following (7.12), 

(7.14) M(t ) :=2t C s (N)e 6cN2 

with v — 7 — e — e'. In what follows, unless stated differently, M will mean M (to). 

First we will bound the linear term in the Q n in (7.3). The sum over ni2 and 
m3 is bounded similarly to the sums in the inhomogeneous term and the sum over 
mi is bounded using the bootstrap assumption (7.13), resulting in 

(7.15) 

N-l 



< 



const 



-N'e 



2-cN 



1 Qmi 1 






r (fe) 




°m 2 




°m 3 



m 2 m 3 0<fc+Z<2(JV-l) mi 
N-l 

<M ^ E < g'+fc+lgc(fc 2 +' 2 ) e c'(fc+0 |ymim 2 m 3 | |x mi | g - (~i-r,'-s') ( | z m2 | + \x m3 \ ) 

m 2 m 3 0<fe+(<2(jV-l) mi 

for | (7 — e') < i/ and ^ < 7 — 77'. This is similar to the sum in equation (7.9) and 
gives a result with a similar dependence on \x n \, 

N-l 



(7.16) s 2 < E /3' +fe+1 Ei^ mim2m3 il^mil 

m 2 m 3 0<fe+i<2(JV-l) "H 

< C { (]V)e 2cAr2 Me"( 7 " £ - E ')l x "l. 

with the same probability as in (7.11). Therefore the linear term in Q n is bounded 
by the probabilistic bound 
(7.17) 

/ ol+k+l ™imjn»3 1 

p J Zjraim 2 m 3 ;!l' J I K n 



> /?C« (iV) e 4cN Me 



-(7- 





c {l) 




r (fe) 




|Qmi | 






c m3 




e—e' — rj' 







with C5 = O (^ 2 )- A similar bound is found for the third term on RHS of (7.3). 
The fourth sum of equation (7.3) 



(7.18) 

is bounded by 
/ 

(7.19) Pr 



S 4 = ]T 

mim 2 m 3 ;fc 



|Qmi I IQ* 



'm 3 



(fc) 



< 



const 



-N A e 



2-cN 
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with C s = (N). The last term in (7.3) is 

(7.20) S a = f3 2N+1 \Vr m2m3 \\Q* mi \\Q^\\Qm 3 \. 

ra\va-2m 3 

and it is bounded by 

(7.21) S 5 < ^AT+lj^g-^-e-e')!^ 

To summarize, substitution of the equality (7.13) in (7.3) to find a time to f° r which 
assumption (7.6) is valid results in the following inequality which is satisfied with 
the probability that is the sum of the probabilities given by the RHS of (7.12), 
(7.17), (7.19), 
(7.22) 

2i Si < to-Sr+toCs (N) (e 6cAf2 + 2[3e AcN2 M + 4(3 N+1 e 2cN2 M 2 + 8/3 2JV+1 M 3 ) e~b~ 

with C s = (N 3 ). Multiplying by e ^~ £ ~ E ') M both sides of the inequality (7.22) 

and taking the infimum with respect to n, gives 

(7.23) 

inf (5 ie (T- £ - £ ')l-"l) < 2C S (N) (/3e 4cjy2 M + 2(3 N+1 e 2cN2 M 2 + 4/3 2Ar+1 M 3 ) . 

Setting A N := e 2cAr2 , C := inf„ ^5ie( 7 - £ - <r ')l 2: "l^ and using the definition of M 
(see (7.14)) 

(7.24) C < 2(3t Q ■ A%C S + Ap N+1 t 2 Q A 7 N C 2 s + 8/3 2N+1 t 3 A 9 



N- 



For sufficiently small (jto the first term on the RHS of (7.24) is dominant and 
therefore, 

(7 - 25) '» - vkc, 

Note, that C > 0, since it is an infimum of a sum of positive quantities, however 
we do not calculate it explicitly in this paper, nevertheless it is likely to be of the 
order of C$A 3 N . This proves that, 

\Q n (t)\ < M (t) ■ e-^- £ - e >" 1 =2t-C s e 6cN2 e~^~ e ~ e '>^ 

for times t < to , where M (t) is given by extending the definition (7.14) by replacing 
t by t. 

Theorem 14. For t = O with c = 2, and assuming that Conjecture 3 holds 

(7.26) Pr (\Q n \ > 2t ■ C s e ecN2 e~^ £ ~ e ')^) < c ^l N ^ e -^ 



The bound on the probability of statement (7.26) is calculated by summing the 
RHS of ((7.12), (7.17), (7.19)). The contribution of the remainder term is 

(7.27) \P N Qn\ < COnst ■ e 6c^+iVln/3 + l„ te -( 7 - £ -s')KI _ 

Note that for a given t and (3 there is an optimum TV for which the remainder is 
minimal. Additionally, for any fixed time and order N, lim^o l/^QnJ / P 14 ^ 1 = 0, 
which shows that the series is in fact an asymptotic one [63]. 
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8. Summary 

In this paper a perturbation expansion in powers of (3 was developed (Sections 
2, 3) for the solution of the NLSE with a random potential. It required the removal 
of the secular terms for this problem. To best of our knowledge it is the first time 
it was done for a multivariate problem. The quality of the expansion to the second 
order was tested in Section 6. In Section 4 it was shown that the number of terms 
grows exponentially with the order. In Section 5 a probabilistic bound on the 
general term (5.29) was derived. It relies on the Conjecture 3. The resulting bound 
was tested numerically. Finally, a bound on the remainder was obtained for a finite 
time, showing that the series is asymptotic. For time shorter than to which is given 
by (7.25) there is a front x (t) oc In t such that for x n > x (t) both the remainder, 
(3 N Q n (t) and c„ (£) are exponentially small. 

The work leaves several open problems that should be subject of further research: 

(1) Turning the perturbation theory developed in the present work into a prac- 
tical method for solution of the NLSE and similar nonlinear differential 
equations. The control on the error should be obtained using the methods 
presented in Sec. 7. 

(2) Can the front x (t) oc In t be found for arbitrarily long times ? 

(3) The asymptotic nature of the series. Is it just an asymptotic series or a 
convergent one ? 

(4) If the series is asymptotic can it be resummed ? 

(5) The e ck +c k in (5.29) results from the repeated use of the Holder inequality 
and a very generous bound on the probability distributions. An effort 
should be made to improve it. 

(6) There are various properties of the Anderson model that have been used 
here. Some of them were tested numerically. It would be of great value 
if Conjecture 3, Corollary 4 and Conjecture 5 were rigorously established, 
even at the limit of a strong disorder. In the present work we rely only 
on Corollary 4 (that was tested numerically, see Fig. 5.1). The rigorous 
proof of the unimodality of 7 (E) for the uniform distribution of the random 
potentials, e x , may also be useful. 

(7) It would be very useful if Conjecture 6 could be rigorously obtained. 
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Complex Systems in Dresden in 2007 and while YK visited the department of Math- 
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Theoretical Physics at the Technion in 2007. The visits were supported in part by 
these Institutes. 



PERTURBATION THEORY FOR THE NONLINEAR SCHRODINGER EQUATION. 



27 




Figure 8.1. The Lyapunov exponent as a function of energy for 
the uniform distribution defined in (5.12). The solid (blue) line 
is for J = 0.25, the dashed (green) line is for J = 0.5 and the 
dot-dashed line (red) is for J = 1. 

Appendix 

The average Lyapunov exponent was calculated numerically using the transfer 
matrix technique for a uniform distribution defined in (5.12). The results are pre- 
sented in Fig. 8.1. It can be seen that 7 (E) is unimodal. 
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